SPI OBSERVATIONS OF POSITRON ANNIHILATION RADIATION FROM THE 4TH 
GALACTIC QUADRANT: SKY DISTRIBUTION 
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ABSTRACT 



During its first year in orbit the INTEGRAL ob- 
servatory performed deep exposures of the Galac- 
tic Center region and scanning observations of the 
Galactic plane. We report on the status of our analy- 
sis of the positron annihilation radiation from the 4*'' 
Galactic quadrant with the spectrometer SPI, focus- 
ing on the sky distribution of the 511 keV line emis- 
sion. The analysis methods are described; current 
constraints and limits on the Galactic bulge emission 
and the bulge-to-disk ratio are presented. 

Key words: gamma-ray observations; positron anni- 
hilation line; Galactic Center region. 



1. INTRODUCTION 



The cosmic positron annihilation radiation was first 
detected through its hallmark, a gamma-ray line at 
511 keV, in balloon observations of the Galactic Cen- 
ter (GC) region in the 1970s and has been the fo- 
cus of intense scrutiny by a large number of bal- 
loon and sat ellite borne experiments ever since (e.g. 
iHarri j [l997j) . Despite this tremendous effort, the 
origin of the positrons is still far from being un- 
derstood. A large variety of positron sources and 
production mechanisms have been pr oposed over the 
years (e.g. IChan fc Lingenfelteilll993l) . Among the 
more promising source candidates are the nucleosyn- 
thesis products fr om supernovae. M ore recently, hy- 
pernovae/GRBs IjCasse et al."2004') and the annihi- 
lation of dark matter (e.g. Bochm ct al. 2003) have 
been revisited. Another intriguing candidate source 
of positrons comprises compact objects. 
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Investigations of the sky distribution of the annihi- 
lation radiation promise to provide clues for identi- 
fying the source(s) of positrons in our Galaxy, de- 
spite the fact that positrons may travel from their 
source before annihilating. The annihilation radia- 
tion was mapped for the first time with the OSSE 
instrument on board the C ompton Gamma- Ray Ob- 
servatory Ce.g. lPurcell et alJl997HMilne et alJ200l|) . 
At least two extended spatial components, a (pos- 
sibly dominant) extended bulge and a disk, are re- 
quired to account for the observed annihilation emis- 
sion. Early OSSE maps showed tantalizing evidence 
for a third emission component (the so-called positive 
latitude enh ancement, PLE) loc ated about 8° north 
of the GC ()Purcellet alJlT997l) : however, the fina l 
analysis did not corroborate this ("Mil ne et alJl200ll) 
and its (astrophy sical and/or instrmnental) origin re- 
mains uncertain l|von Ballmoos et al.ll200,'^ . 

Modelling of the Galactic distribution of the 511 keV 
line emission from OSSE data suggests FWHM val- 
ues of about 4° — 6° when using Gaussians to describe 
the bulge compo nent. No significant off set from the 
GC w as found by Purcell et al. (1997) or 'Milne ct alj 
l|2001f) . However. Tueller et ah (|1996) . using OSSE 
observations with the long axis of the collimator per- 
pendicular to the Galactic plane to produce a maxi- 
mum sensitivity scan in longitude, derived a signifi- 
cant ('^ A(j) offset of the 511 keV line emission. An- 
alyzing d ata from the wide field-of- view TGRS in- 
strument, iHarris et al.l l)1998D obtained a rather ex- 
tended distribution (FWHM ^ 25°) with no signifi- 
cant offset from the GC. Before SPI, OSSE was the 
only instrument that allowed the bulge-to-disk ratio 
(hereafter B/D) to be constrained, albeit poorly, to 
a range of 0.2 - 3.3. This parameter is important for 
characterizing the Galactic distribution of the emis- 
sion and - potentially - of positron sources. A major 
uncertainty in the determination of B/D is the pos- 
sible existence of an extended (halo) component to 



the bu lge, which w ould lead to a large ratio. Using 
OSSE. lMilne et alJ l|2000|) could determine the total 
511 keV line flux from their model Galaxy relatively 
well (unlike B/D): (2.1 ~ 3.1) x IQ-^ ph cm'^ s'^. 

Early SPI results on the 511 keV line em ission from 
the GC re gion have been reported by iJean et all 
l|200.1t) and iKnodlseder et all i|200.'^ . The Galactic 
distribution of the emission could be described with 
a Gaussian bulge component of about 9° FWHM and 
a 2a uncertainty range covering 6° — 18°. The cen- 
troid of the Gaussian was found to be consistent with 
the GC at the 2.1a level. No evidence for a Galac- 
tic disk component was found. However, flux limits 
were in agreement with OSSE measurements. 

In this paper we report results on the sky distribu- 
tion of the 511 keV positron annihilation line radi- 
ation from the GC region after one year of obser- 
vations with the spectrometer SPI on board ESA's 
INTEGRAL observatory. Spectroscopic analysis of 
th e 511 keV line , using the same data, is presented 
bv iLoniou et al.l l)2004aj) . The results from both of 
these studi es are prese i ited in the context of previous 
findings bv lJean et al.l 1120041): their astrophysical im - 
plications are discussed bv iG^iessoum et all (j2C 



2. DATA ANALYSIS 



The data used for the analysis of the sky distribu- 
tion of the 511 kcV line emission with the spec- 
trometer SPI l^drcnnc et al. 2003) consist of the 
spring and autumn Galactic Center Deep Exposures 
of 2003 (hereafter GCDE), supplemented by Galactic 
Plane Scan (GPS) observations of the central Galaxy. 
Both the GCDE and the GPS observations are part 
of the so-called Core Program observations of the 
INTEGRAL observatory, which are proprietary to 
the INTEGRAL Science Working Team (ISW T) an d 
the instrument teams for one year ( Winkler! l2001(l . 
The data comprise the INTEGRAL orbital revolu- 
tions 47-66 (March 3 - April 28) and 97-123 (Au- 
gust 2 - October 16), combining 1266 and 1068 point- 
ings with a total exposure time of about 1.8 and 
2.0 X 10^ s, respectively. The GCDE consists of rect- 
angular pointing grids covering Galactic longitudes 
|/| < 30° and Galactic latitudes |6| < 20°; the GPS 
consists of poi ntings along th e Galactic plane within 
|6| < 6.4° fsee lWinkledl20(nL for details). 

Due to data sharing agreements within the ISWT, 
the results presented here are limited to the 4*'' 
Galactic quadrant {I = 270° — 360°). However, in ac- 
cordance with the aforementioned agreements, data 
from pointings in the entire GCDE region (as de- 
fined above) have been included in the analysis. The 
resulting exposure to the sky is depicted in Fig. ^ 

The data were prepared following the procedures de- 
scribe d in Ijea n et al. ( 20(^ and [Knodlscdcr c t al.l 
l|2003D . SPI single detector event s have been gain 
corrected (see ILoniou et al]l2004b() and binned into 
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Figure 1. The sky exposure resulting from the Core 
Program observations used in the analysis 



event spectra for each detector and instrument point- 
ing with 0.5 keV bins, leading to a 3-dimensional 
data space. In this data space, the instrumental 
background in the 485-550 keV band was modelled 
by two components: the instrumental 511 keV back- 
ground line and the underlyin g continuum back - 
gromid, as described in detail in 'Jean et al.' (zQQ^)) 
tTccgardcn ct al. (200|), andpcan ct al. (2004). 

The continuum background in the line analysis in- 
terval 508.5-513.5 keV was determined based on the 
continuum background level at adjacent higher ener- 
gies. The modelling of the instrumental line compo- 
nent has b een significantly im pro ved since the first 
analys es bv l.Iean et al.l ^^2003^ and IKnodlseder et al.1 
l|2003t) . We employed a multi-component model for 
the instrumental 511 keV line, which reflects to first 
order the multitude of prompt processes and delayed 
(radioactive) decays that give rise to this background 
line. We used e.g. the SPI background line iden- 
tifications l|Weidenspointner e t al. 200^ and instru- 
mental background simulations employing the MG- 
GPOD Monte Carlo suite (Weidcnspointncr ct al] 
l2004albf) to identify the main contributors to the 
511 keV line. The prompt components were assumed 
to scale with the rate of saturated Ge detector events 
(GeDsat). Delayed components due to radioactive 
decays were modelled by calculating a time history 
of their activity, assuming that their production rate 
scales with GeDsat. 



3. RESULTS 



Imaging analyses of t he 511 keV line emission with 
Richardson-Lucy (see IKnodlseder et al.ll2003[ and 
references therein) and Bayesian l|Allain et al.ll200^ 
algorithms clearly show extended emission from the 
Galactic bulge with a radial profile that is approxi- 
mately Gaussian. The sky maps do not provide any 
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Figure 2. The best fit model for the 511 keV line 
emission from the GC region. A Gaussian, repre- 
senting the Galactic bulge, is sufficient to explain the 
data. Details are given in the text 



evidence for emission from point sources, a PLE com- 
ponent, or (at this stage) the Galactic disk. 

A more quantitative approach for studying the 
Galactic distribution of the emission is model fit- 
ting. As a first step we have modelled the bulge 
emission by a spherical distribution with a Gaus- 
sian radial profile, located at the GC. We obtain a 



best fit FWHM of 8° 
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with the quoted uncer- 



tainties defining the statistical 2cr confidence lim- 
its. The corresponding bulge flux is (0.961q'^4) x 
10^'^ ph cm^^ s^^, with the uncertainty being domi- 
nated by the uncertainty of the width of the Gaussian 
intensity distribution. Spectroscopy of 511 keV line 
emission from the bulge resulted in a best fit energy 
of (511.02to'o9) keV and an intrinsic fine width of 
(2.67+J];|!^) keV FWHM (see lLoniou eTani2n04aD . 

If we assume an ellipsoidal distribution with a Gaus- 
sian radial profile we improve the fits marginally. A 
bulge distribution that is more extended in longitude 
than latitude is slightly, but not yet significantly, 
favoured. This elongation might be an indication 
of emission from the Galactic disk. 



per limits). The DIRBE 35 /im and 240 /im maps, 
which trace old and young stellar populations, yield 
flux limits of 1.4 and 0.9 xlO~^ ph cm^^ s^^ when 
integrated over the Galaxy. The motivation for chos- 
ing these two models was to span plausible positron 
source distributions. 

We tried to account for the 511 keV line emission 
with point sources only using both model fitting and 
SPIROS (Skinner & ConneU 2003) imaging analysis. 
A single point source can be excluded with a very 
high confidence level. However, formally we cannot 
exclude the possibility that the emission originates 
in a population of at least 4 point sources. We also 
investigated how much a single point source could 
contribute in addition to our best (diffuse) descrip- 
tion of the bulge emission. For a point source at the 
position of Sgr A*, we can already set a 2cr upper 
limit of about 40% of the total fiux. 

Analyses of the putative PLE component are compli- 
cated by the fact that the various OSSE studies did 
no t converge on a unique pa rametrization (see table 
in Ivon Ballmoos et al.l 12003(1 . In the following, we 
restric t ourselves to a PLE model of Purcell et alj 
l|1997t ). a Gaussian at I = -1.8° and b = 11.6°, 
with a FWHM of 16.4°. Fitting this PLE model to- 
gether with a Gaussian model for the Galactic bulge 
has no significant effect on the results for the lat- 
ter. No significant flux is attributed to the PLE; a 
conservative 2a upper limit on the PLE fiux (taking 
into account the 2cr uncertainty in the bulge size) is 
1.6 X lO^"* ph cm^^ s^^. It should be noted that 
a PLE of the strength suggested by iPurcell et al.l 
l|1997t) would have already been noticeable in the sky 
maps. 

We have begun to investigate the positronium (Ps) 
continuum emission. First results in the 409-505 keV 
band indicate that the sky distribution of the to- 
tal emission is consistent with what we find for the 
511 keV line. If we use the lHarris et alJ l|1990^) result 
to fix the relative intensities of Galactic diffuse and 
Ps continuum emission, we can infer a Ps continuum 
flux from the measured total flux that is consistent 
with previous findings. However, because of system- 
atic uncertainties in the subtraction of Galactic con- 
tinuum emission we refrain from quoting a value for 
the Ps continuum fiux and the Ps fraction until a 
forthcoming publication. 



If we do not constrain the spherical Gaussian model 
(FWHM 8°) to be located at the GC, we find a best 
fit position of ; = -1.0° ± 0.7° and b = 0.3° ± 0.7° 
(Icr error bars). The best fit bulge model is depicted 
in Fig. [21 The significance of the displacement from 
the GC is marginal at 1.6a. This displacement of the 
Gaussian slightly reduces the 2a range for FWHM, 
and slightly increases the bulge flux. 

Although model fits yield hints of 511 keV line emis- 
sion from the Galactic disk, such a component has 
not yet been significantly detected. The fiux at- 
tributed to the disk depends strongly on the assumed 
spatial distribution (in the following, we quote 2a up- 



4. DISCUSSION 



With the available data, the sky distribution of the 
511 keV positron annihilation line emission from the 
GC region as measured by SPI can be adequately 
described by a spherical distribution with Gaussian 
radial profile and a FWHM of S°t-2o , the quoted un- 
certainties being statistical 2a confidence limits. The 
bulge extent as determined with SPI is consistent 
with the best fit OSSE results (4° - 6°) at the 2a 
confidence level. We cannot draw any firm conclu- 
sions regarding a possible discrepancy until we have 



achieved a positive detection of the Galactic disk 
with SPI, which might reduce the size of the bulge 
component. 

We detect no significant deviation of the bulge dis- 
tribution from spherical symmetry. The data hint 
at a possible elongation in longitude, which could 
reflect faint emission from the Galactic disk. Allow- 
ing the location of the Gaussian bulge model to vary 
freely suggests a slight offset from the GC that is 
only marginally significant. However, it is interest- 
ing to note that our best fit position coincides al- 
most exactly with the centroid of the 511 keV line 
emission found by iTueller et"an l|199 (JD using OSSE 
data. If the existence of an offset of the centroid of 
the 511 keV line emission can be corroborated in fu- 
ture analyses including more data, this would set an 
important constraint on emission models - in partic- 
ular on models that attribute the bulk of the bulge 
511 keV line emission to the annihilation of light dark 
matter particles, whose distribution should be ex- 
actly aligned with the dynamic center of the Galaxy. 
It is also interesting, but perhaps coincidental, that 
the centroid of our best-fit Gaussian bulge distribu- 
tion lies less than la away from the hard X-ray source 
1E1740. 2-2942, which has been reported to produce 
transient red-shifted a nnihilation line radiation (see 
review bv lHarriall997l and references therein). 

Assuming that the centroid of the 511 keV line emis- 
sion is the GC, we find a flux from the bulge of 
(0.96t!^;i4) X 10-3 ph cm-2 s-\ with the uncer- 
tainty being dominated by the uncertainty of the 
width of the Gaussian intensity distribution. This 
value is in good agreement with previous results ob- 
tained by high-r esolution instruments (see table in 
iJean et all |2003() , which had flelds-of-view smaller 
than or comparable to that of SPI, but falls short 
of the total 511 keV line fluxes obtaine d with OSSE 
fe.g. Milne et al. 20001 and SMM (e.g. Karris et all 
119901) bv a factor of 2-3. SMM had a very large 
field-of-view (~ 130° FWHM) and therefore was sen- 
sitive to very extended emission beyond the Galac- 
tic bulge. In case of OSSE, the fiux was sepa- 
rated into bulge and disk components, which had 
the values (0.5 — 2.4) x 10^"^ ph cm^^ s~^ and 
(0.8 — 2.6) X 10"'^ ph cm^^ s"^, respectively; these re- 
sults strongly dep ended on the assumed shape of the 
bulge component (Mi lne et al.l l2000'l . Our results for 
the bulge flux, and our 2a upper limits on the disk 
flux, (0.9 — 2.2) X lO^'^ ph cm^-^ s~^ , are consistent 
with the OSSE measurements. 

Our upper limits on the fiux from the Galactic disk 
derived from model fits using the DIRBE 35 ^m 
and 240 /im maps can be converted into lower lim- 
its on B/D, which are 0.5 and 0.8, respectively. 
Again, these limits are compatible with OSSE mea- 
surements (0.2-3.3). The SPI limits in B/D already 
have interesting implications. Based on their own 
low B/D, both DIRBE maps can be excluded as sole 
source distribution of Galactic positrons. Similarly, 
B/D for Galactic Type la SNe (from the distribu- 
tions by iDawson fc JohnsonI Il994t iMatteucci et al.l 
l2003|) is below the SPI limits. However, the observed 



distribution of low-mass X- ray binaries is consistent 
with our results for B/D l|Pranzosl l200^ . If light 
dark matter particles provide an important source of 
positrons, these would have a high B/D value. 

We do not find any evidence for a PLE component. 
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